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We analyze the quantum statistical treatment of bound states in Hydrogen considered as a system of electrons 
and protons. Within this physical picture we calculate isotherms of pressure for Hydrogen in a broad density 
region and compare to some results from the chemical picture. First we resume in detail the two transitions 
along isotherms : 

(i) the formation of bound states occurring by increasing the density from low to moderate values, 

(ii) the destruction of bound states in the high density region, modelled here by Pauli-Fock effects. 
Avoiding chemical models we will show, why bound states according to a discrete part of the spectra occur 
only in a valley in the T-p plane. First we study virial expansions in the canonical ensemble and then in the 
grand canonical ensemble. We show that in fugacity representations the population of bound states saturates at 
higher density and that a combination of both representations provides quickly converging equations of state. 
In the case of degenerate systems we calculated first the density-dependent energy levels, and find the pressure 
in Hartree-Fock-Wigner approximation showing the prominent role of Pauli blocking and Fock effects in the 
selfenergy. 
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1 Introduction 

In 1911, exactly hundred years ago, Rutherford invented a new model of matter, which is the basis for the science 
deaUng with strongly coupled Coulomb systems HI. His model was based on the interpretation of existing 
experiments on the scattering of electrons on matter. In May 1911, Rutherford came forth with a model for 
the structure of atoms which provided a first understanding why electron scattering can so deeply penetrate into 
the interior of atoms, so far unexpected experimental results. Rutherford explained his scattering results as the 
passage of a high speed electron through an atom having a positive central charge +Ne, and surrounded by 
a compensating charge of N electrons lUl. In Rutherford's model the atom is nearly empty, it is made up of 
a central charge and electrons and the glue keeping the system together are the Coulomb forces. Soon after, 
in 1913, a first quantum-mechanical treatment of the model was given by Bohr and a statistical treatment by 
Herzfeld. According to the Bohr-Herzfeld theory the energies and the corresponding partition function are 



s=l 



We denote the main quantum number by s in order not to be mixed up with the density n. According to the 
Bohr model there are infinitely many levels close to the series limit s ^ co the terms increase as and the 
sum diverges. This is the result for hydrogen, but the problem remains for all elements. The problem of the 
divergency of the atomic partition function took a long time before the mathematical background for a serious 
treatment of this problem was available. However urgent need to solve the ionization problems in astrophysics 
and plasma physics forced fast solutions. For the calculations the problem was simplified and the partition 
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function was estimated by the first term of the sum, or by two or three terms. Ignoring the open problems with 
the partition function, Eggert and Saha attacked the urgent problems of astrophysics, considering only the ground 
state contribution. 

An idea developed later was to cut the sum at the minimal terms \Es\ ~ ksT. This way of solving the 
divergence problem is due to the classical paper of Planck f2] published in 1924. Developing Plancks approach 
Brillouin proposed in 1930 a more smooth way of removing the divergence which led to the formula ||3| 

oo 

<tbpl(T) = A<^MPI/s^) - 1 - {1^1 1 s^)] ■ (2) 

s=l 

Using semi-classical methods, Planck could show that the (infinite) interaction terms below and above the series 
limit compensate each other. This procedure was justified by methods of quantum statistics only in 1960 by 
Vedenov and Larkin ID |5| and subsequently by the present authors in collaboration with D. Kremp by using 
quantum-statistical methods ||6l|7l|8l|ll[l0l[II][l2l[ni[H[l5l. We will show here that the Brillouin-Planck- 
Larkin procedure provides the most natural way to avoid the divergencies for non-degenerate systems. Further 
we study degenerate systems and show that a consequent treatment of the identity of electrons in the Hartree-Fock 
approximation including bound states is a key for the study of hydrogenic bound states. The present approach 
is based on the method of Greens functions as developed in ifTTl [121 [131 [T4l [TSll and recently surveyed in a book 
dedicated to Nevill Mott fTS). The present investigation is concentrated on the so-called physical picture which 
stays within the Rutherford picture and avoids the introduction of chemical species like atoms, molecules etc. In 
some of the earlier work the transition to a chemical picture was made following the principle of equivalence that 
bound states are to be treated on the same footing as free particles |[T4l[r7ll . This approach was quite successful 
in the description of partial ionization ATI [TSl [191 [20l [211 [22l . On the other hand the chemical picture may lead 
to principal difficulties e.g. if at high densities the minimum of the free energy is not well defined 1231 [24l . This 
is the reason why we stay here within the physical picture trying to sum up the important higher order terms, this 
way following the line developed in several recent works 1251 [26l [27l [28l . 



2 Density and fugacity expansions for nondegenerate systems including non- 
linear ring and Saha-type contributions 

Any successful description of plasmas has to go beyond the Saha theory which is something like the zeroth order 
approximation for plasmas. The derivation of Saha-type approximations from quantum statistics was first studied 
by Planck and Brillouin ||2l [31 and is since then a topic of permanent interest due to the importance of Saha- 
type equations for experimental work and for many technological applications |6|. We use here a systematic 
quantum statistical approach to the pressure on the basis of the Greens function representations of the pressure 

1 d\ f 

p(/3,Aie,/^.) =Pid- ^ — J dldlV{n)G2{l,l,l++,l+ ■.h=ft) (3) 

where 1 = {pi, cri} denotes momentum and spin variables. The Green's function representation works in the 
grand canonical ensemble which provides us a series in the fugacities (s^ - spin of particles, a — e,i) 



Ag y/2'!TmakBT 

Note that is the excess part of the chemical potential (the part beyond the ideal Boltzmann term) and that 
therefore Za ria for small densities. A representation by a density series may be obtain by inversion, i.e. by 
excluding the fugacities. We summarize the existing results about the pressure of a hydrogenic quantum plasma 
as a function of the density in a rather special form which is asymtotically valid in the region of low and moderate 
temperatures T < I /ksT. Using for this region the asymptotic expressions of the exact results, the pressure may 
be written as (k^ = gTr/Sne^ lTTl[T3l 

/3p(/3, n) ^ n[l + ^nA^ + 0{n^)] + n - ^[l - ^(^A) + ^(^A)^ + ...] 

~n^-i[l - fie^n] + 0{n^ log n) . (4) 
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We introduced here the mean thermal wave length 

^ = A,e = „ , ^ = 7-^; M= : ; 7 = AgCTBPL r = 87rV7rA,;eCTBPL r 

2fikB-L 2-y/7r me + rrii 

where A is an effective quantum length, being in some correspondence to the Debye-Hiickel parameter in the 
classical theory. We note as a special property of the given temperature regime that he interaction terms in the 
pressure depend only on reduced mass fi of the electron-proton interaction and not on the other reduced masses. 
An equivalent form of eq.© which shows better the general structure is the following 

/3p(/3, n) = n[l + -^nA^ + 0(n^)] +n- nj[{l - Pe^ n + 0{k^)) + ...] 
Sv 2 

-^[1 ^ ^7 + 0(7^)][1 - ^(.A) + + 0(.3)] + 0(n3/^) . (5) 

The expression for the pressure of hydrogen in the limit of small densities, given in the formulae eqs. (|4| and (|5]) 
has been drawn in Fig. [T] This is a first step to a Saha-type description, since the pressure related to the ideal 
pressure is decreasing with the density with a factor proportional to the distribution function, however we see that 
our formula fails at higher densities going even to negative values. Looking at the curves in Fig. [T]we see that the 
screening effects and the contributions from the Brillouin-Planck-Larkin function both tend to lower the pressure. 
This lowering is too large at increasing densities, so we have to look for contributions limiting this growth. In 
order to proceed in this direction, we investigate first the structure of eq. |5] aiming to complete the unfinished 
series with respect to 7 and kA. The structure of eq. (|5]l suggests that the expressions in the parenthesis may 



Fig. 1 Left panel: Pressure related to the total classical pressure 2nkBT at constant temperatures T = 15788 K and 
T — 23682 K in dependence on the density (log-scale). The result of density expansions up to 2nd order virial terms (red and 
green) are compared to corresponding fugacity expansions. Right panel: Pressure related to the full ideal pressure (including 
the Fermi pressure) calculated from an extended density representation including the ring sum ^(x) and the Saha function 
A(x) for the temperatures T = 15000 K, 20000 K, 50000 K in dependence on the density (in log-scale 15 - 25). 

be only the first terms of some infinite series with respect to nA^, k and 7. The first series is evidently nothing 
else than the ideal electron pressure, which may easily be extended to an infinite series using the ideal Fermi 
pressure. In order to find the higher orders in the series in 7 and in kA we will use a more complicated procedure 
which is based on fugacity expansions lfTn[r7l|29l According to our general philosophy, the strong decay of the 
pressure observed in second order density expansions (Fig. [T]) is due to the fact that important higher order terms 
corresponding to some infinite series were omitted in the derivation based on the canonical ensemble. Alastuey 
et al. II25II26II succeeded in deriving a Saha type expression in summing up an infinite series of terms. Here we 
follow a similar idea but go an alternative way through fugacity expansions. As shown much earlier, fugacity-like 
expansions are equivalent to mass-action laws. This was used in plasma theory by Bartsch, Ebeling and others 
IfTTl [TTl |28l |29l and worked out in large detail by Rogers et al. 1001 . We mention that Rogers et al. succeeded to 
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show that extended fugacity expansions allow an excellent description of measurements for the oscillation modes 
of the sun 1311 . 

A systematical quantum statistical approach to the pressure including the low as well as the high pressure region 
may be given on the basis of Green's function representations of the pressure ifTTl [121 [131 [T4l . The Green's 
function representation works in the grand canonical ensemble which provides us a series in the fugacities instead 
of a series in the densities. Note that jj!^^ is the excess part of the chemical potential (the part beyond the ideal 
Boltzmann term) and that therefore Za — > ?ia for small densities. The structure of this series is similar to the 
density series, we have e.g. similar screening effects, however there are some differences, e.g. the fugacity 
series contains more diagrams. The contributions of bound states of an electron and an ion (atoms) are contained 
in the contribution of order ZeZi, the contribution of molecules is contained in the term of order z^zf. General 
expressions for screened fugacity series were written down first by Montroll and Ward, Vedenov and Larkin and 
explicit calculations were done by De Witt and Larkin ||4||5][32l. Semiclassical fugacity series were given by 
several workers OTl [29l . We are interested here mostly in the bound state effect and may simplify by using 
approximations allowed for the region T < 157000 K were the plasma is nondegenerate and where the bound 
states are dominant. As we have seen in the density representation, the differences of the masses of electrons and 
ions (protons) do not play a role in the region T <^ I kg and the plasma behaves asymptotically like a system 
with equal masses and equal relative De Broglie radii. In the lowest approximation corresponding to eq. (j?]) we 
find for low temperature H-plasmas the following results flTl [131 [33l 



Ze+ Zi + 



—^f{ng\) + STTZgZjA^ exp[l 



^e^Kg]<jBMT) + 0{z 



5/2N 



(6) 



where the grand-canonical screening quantity is defined by — 47r/3(ze + Zije^ and the quantum-statistical ring 
function is a transcendental function having the series 



1 



1 



16 10 

The densities may be expressed by fugacities 



d 



-r^ = z, + -l — {KgXf{Kg\)) + A7TZeZ,[l + j3e^Kg + • • • ]A3(Tbpl(T) ■ 



16 dX 



(7) 



(8) 



< d 



-^=Ze+ T^l^iKgXfiKgX)) + ATrZeZ,[l + /3e^Kg + • • • ]A^CrBPL(T) 



dz, 



16 dX 



(9) 



We will show by iterations that going from the fugacity variable to the density this result occurs to be equivalent 
to eq. (jUi. In order to proceed with this complicated system of equations we go to the special case of non- 
degenerate hydrogen at lower temperatures T < I /ks and satisfying j^A^ < 1. As shown above, in this region 
the differences of the electron ion masses does not matter and we my assume z^ = Zi = z. In order to represent 
the pressure by densities we neglect now in a first step the bound state contributions assuming crBPL(r) = 0. 
Expressing in this approximation first the fugacities by densities we find 

„2 "1 



n exp 



G{kX) 



G{x) = ^ 

X 



1 



exp(^) 



1 - $(-) 
2 . 



= 1 - 



1 2 

6" 



, (10) 



where <I>(a;) is the error function. Then excluding the fugacities from the series step by step we arrive at the 
representation 



f3p = 2n- -^Hi^^)'^ 



(11) 



Here (j>{x) is another transcendent function, related to f{x) and f'{x), the so-called pressure ring function which 
is an analogue of the Debye-Hiickel ring function 141 [iTI . A full representation of this transcendent function is 
given by the infinite series 



3 



fc=i 



2fefc! 



It. 



k=l 



(fc + l)(fc + 3)a; ' 
2fc(2fc- 1)!! 



2k 



(12) 
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In the calculations we used more convenient Pade formulae. The nonlinear functions 4>{x) and G{x) converge to 
zero at large arguments what is important for the high-density asymptotics of the pressure. In Fig. [T]the result of 
fugacity expansions up to 2nd order virial terms for T — 15788 K and T = 23682 K (red and green curves) are 
compared with corresponding density expansions. We see that the behaviour of truncated fugacity expansions 
seems to be much better than that of truncated density series. However this is not true in all cases and may lead 
to serious problems. We may show this on the example of the ideal Fermi gas. Truncating here the series after 
quadratic terms the representation of the pressure leads for nA^ > \f2 to an imaginary pressure. For this reason 
we hold the opinion that the a controlled summing up of terms in the density series may be more reliable. The 
results we obtained so far may be summarized as follows. By means of grand-canonical methods we succeeded 
already to obtain one infinite series in the screening parameter k instead of the first linear and quadratic terms 
only we found in the density series, (see eq.©). This summing up a series in k which leads to a saturating 
function what improves very much the behavior at larger densities. 

In our canonical representation eq. (|5]) appears also another series expansion in the density parameter 7. In order 
sum up this series in 7, we proceed in a similar way. We study at first the bound state contributions, neglecting the 
contributions coming from the ring terms and from the terms due to degeneracy. This yields the simple quadratic 
equation 

n = z + 87rz2AVBPL(r). (13) 

By the way, this formula shows already that the fugacities of the electrons and the protons should be rather small 
in the bound state region where the partition function is large. The quadratic equation for the fugacities can be 
easily solved. This way we find after elimination of the fugacity z\ 

z°==nA(7); A{x)^^yYTVx-\\- -f ^ SirnX^ aBPL{T) . (14) 

2x 

Introducing the zeroth step of iteration 0° into the pressure we find 

f3p° ^n{l + A{j)). (15) 

This representation includes a nonlinear function A{x) which saturates at large densities and is closely related to 
the result of the Saha theory, except that in our version the BPL partition function appears. The new representation 
contains all terms in the density up to as well as several higher order terms in the density contained in the 
nonlinear Saha function ^(7). We will show that this kind of function allows us to reproduce an ideal Saha-type 
behavior. We note the following series expansion and asymptotics of the nonlinear function A{x): 

A{x) = l-x + 2x^ ; A{x) ~ ^ + • • • . (16) 

The saturating behavior of the nonlinear function A{x) in comparison to its linear approximation has the conse- 
quence that in the region of large partition functions the fugacities disappear as 

Z° = =. (17) 

This is important for the understanding, why the fugacity series has in the region of bound states a better conver- 
gence as the density series. The better convergence is due to the fact that the fugacities disappear in the region 
where the bound state contribution is large. The fugacity series contains more terms than the density series, so we 
may expect that some of the difficulties connected with density expansions, as the strong decrease of the pressure 
with increasing density shown in Fig. [T]may be avoided. Indeed, a representation of the curve corresponding to 
eq. (|6]l shows a more reasonable behavior with increasing density which is much closer to the Saha-type behavior. 
As the curves shown in Fig. [T] demonstrate, the pressure according to the fugacity expansion goes to saturation. 
This is due to the fact that the fugacity is not fixed, it decreases with increasing value of the partition function and 
this way limits the growth. Evidently the fugacity expansions provide a more correct description of the bound 
state contributions as demonstrated in Fig. [T] In order to draw some first conclusions; 

Density as well as fugacity expansions have both advantages as well as disadvantages: 1) The density expansion 
describes well the screening effects but it fails to cope with the diverging contributions from the screening terms 
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and from the BPL-partition function. 2) The fugacity expansion corresponds to an infinite density series including 
the partition function a. If a is large, then the fugacity goes to zero what guarantees always finite contributions 
to the pressure. This is also true for the screening contributions; any strong increase of contributions damps the 
fugacities. 

We may expect that the best representation is obtained by extended density expansions which contain additional 
contributions corresponding to the important damping terms in the fugacity expansions and provide Saha-like 
terms. This kind of mixed expansions combine the positive features of both expansions avoiding the negative fea- 
tures. In order to demonstrate this we started here from the density expansions and used the fugacity expansions 
for finding the right continuation of the infinite series with respect to the 7- parameter and the k\ parameter. 
This way we obtain in a first step the following expression for the fugacity and the pressure 

z^nA (cxp[-/3e2KG(KA)]7) , (18) 



/3p = n - -^<?^(«:'^) + (cxp[-/3e2KG(K'A)]7) , k! = tiMafl"^, (19) 

where a renormalized value of the screening parameter k! appears. This provides us a closed and relatively simple 
formula for the pressure which however is restricted to the nondegenerate region. However so far our results refer 
to non-degenerate electrons. The good convergence of the new expression including nonlinear functions A(x) 
and ^(x) based on the fugacity series is explained by the fact that the fugacities of the electrons and the protons 
are rather small in the bound state region. The overall behavior of the new representation is much better than 
that of pure density or pure fugacity representations. We may conclude that the most appropriate description of 
Coulombic systems is by density expansions and including infinite sums representing screening and bound state 
effects. We note that the new theory based on eq. (fT9] l is consistent with the Saha-theory and in particular also 
with the Saha-Debye-Hiickel theory ifTTI . We underline that several of the terms beyond v? and are based on 
extrapolations which still need further confirmations. However in the region of low temperature T < //fc^ and 
non-degenerate plasmas our rather simple formulae give a rather good behavior and describe well the transition 
from low density to the valley of bound states. We note however that several physical effects as e.g. plasma phase 
transitions ifTTl [TSl are not yet described by the present approach. Maybe this effect appears only after transition 
to some chemical picture ifTTl ITSi . 

In order to describe as well the transition to full ionization in the degenerate region, additional effects are to be 
taken into account, in particular the symmetry between electrons and protons is lost. The most important new 
effects at high densities are connected of the replacement of the ring terms by Hartree-Fock terms and a density 
dependence and final disappearance of the bound states with increasing density. 



3 Calculation of the pressure at high densities in mean field approximation 

The theory presented in the previous section provides isotherms of the pressure for the nondegenerate region. 
After crossing the line of degeneracy riA^ ~ 1 the theory is no longer valid since degeneracy effects and the 
dissymmetry of the massed are relevant. In this region the formulae given above are just an extrapolation which 
is exact only in the (dominant) ideal Fermi contributions. For very dense, strongly degenerate plasmas, bound 
states do not exist, the plasma behaves like a degenerate nonideal Coulomb gas. Let us assume, we start from 
a high density and decrease the density along an isotherm. The question is then, at which density bound states 
will appear. In other words we have to ask the question at which densities the spectrum will have discrete levels 
corresponding to bound states of two particles, three particles, four particles etc. Turning the question around 
we may ask, what is the density dependence of discrete states. This way, the most important effect to be taken 
into account in the degenerate region nA^ > 1 is the change of the energy levels by effective density-dependent 
energy levels and the question of merging of discrete levels in the continuum. This strongly affects the equation 
of state of hydrogen lElEl- 

In the region of strongly degenerate electrons, the thermodynamics is in good approximation described by the 
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Hartree-Fock approximation, the protons form a kind of a Wigner lattice ITSl . Based on eq. ^ the grand- 
canonical pressure may be approximated by lfm[T3l[T4ll 



p{l3,ze,z,) = p'''+pf''{z,)+pf^'{z,) + z,z,Ald{T) + ... (20) 



(21) 



with 



Ze + Ze r h ZeZiA^Cr(T) , (22) 



d 



Ze 



n, = z, + z,^-^^^-^ + ZeZ,Ala{T) . (23) 

We proceed as in the previous section. Neglecting first the bound state contribution we get for the fugacity of the 
electrons in Hartree-Fock approximation ifTTl [T3l 

g2 

Ze = nexp[-^— /_i/2(a)]; a = (24) 

Here the Ik{x) are the well-known Fermi integrals. We use in the following the high density limit including 
the first temperature corrections calculated with the Sommerfeld method ifTTl [121 [131 [35]| (t = fc^r/Ry). Since 
the correction in the Sommerfeld expansion changes the sign with increasing temperature we used in the 
calculations a more convenient Pade expression with the T^— term in the denominator 



0.3059 n „p , 1.222 



Pf^ [R-y] , non ^2 ; zf^ c^n exp 



irs +0.1216^3 r^' 



-,1/3 



Airn 



■ (25) 



The ionic contribution at high density is more difficult fT3\. However in the limit of high density, the Wigner 
limit, we may give the result in an analytically quite similar form with lowest order T-corrections: 

. 2.3856 (1 + 0.05 r,), wiru i 0.5964 ?i(l + 0.01 r,) 
Zt ~ nexp ; pj Ry ~ . (26) 

In difference to the situation at small densities, the fugacities are at high densities (small r^) rather small in com- 
parison to the densities. For this reason the contribution of the bound states to the pressure is just a perturbation 
and we do not need to sum up an infinite Saha-type series, as in the low density case. Putting now together 
all terms we arrive in the canonical ensemble at the following simple relation for the high density pressure of 
hydrogen (p^ - Fermi pressure of the electrons) 



p(/3, n)c^pl+ pf"" + pY' + z^^'zT'Ki [^M-PEr.) 1 + PEr. 



(27) 



This way we found now for the limit of high density a closed expression for the pressure of hydrogen in the 
physical picture. It remains to find the effective two-particle energy levels which are now density-dependent. 
The two-particle bound state energies in a dense system £"„ may be derived from an effective wave equation, the 
so-called Bethe-Salpeter equation. 

In this paper, the influence of the medium is taken only in mean-field approximation, i.e., by Pauli blocking and 
Fock self energy . The effective wave equation reads in this case (in the following the tilde is omitted) |[T3l[34l[36l 

^^„(p) + ^ F((z)Vn(p + 9) - E^i^^ip) = ^('?) [^"(^^ + - Mp)fe{p + q)] ■ (28) 

Here the l.h.s. represents the standard Schrodinger equation in momentum representation. The terms on the 
r.h.s. containing the electron distribution function fe{p) model effects connected with the Fermi character of 
the electrons. The first of these two contributions we denote as Pauli-blocking and the second one as Fock 
contribution. The Fock term and the Pauli blocking contributions have opposite signs and compensate each 
other partially. The more complete Bethe-Salpeter equation contains also dynamic self-energy and dynamically 
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screened interaction contributions. These corrections were studied for hydrogen and electron-hole systems, see, 

e.g., IE] [Ml [361 [321 EH im Him]. 

The corresponding shifts will be neglected here since we consider only the high density limit. In perturbation 
theory, the sum of Pauli and Fock shifts ITSl [39] l42l can be calculated with the solutions 4>n of the unperturbed 
Schrodinger equation. 

Ai?r = - E Mp)viq)[feip)Mp + q)- feiP + q)Mp)]- (29) 

P,<1 

In the limit of strong degeneracy where the Fermi function can be replaced by the Heaviside step function, the 
integrals can be performed analytically. The result for the shift of the lowest three energy eigenvalues is shown 
in Fig. [2 In the low-density limit we have for the ground state = {10} 

lim A^Pf = 77 y ^0io('?)[<^io(O) - <i>M\ = 327rn' - 207rn' = 12^7i' . (30) 

9 

In analogy we find analytical expressions for the excited s-states. In the low-density limit we have 

lim l^El^ = 48?™' , lim AE^J = IOSttji' . (31) 

n— >0 n— ^0 

As we see, due to compensation effects, the sum of Pauli and Fock shifts is, in linear approximation, less than a 
half of the Pauli shift in the ground state and is also reduced in the excited states. The evaluation of eq. ( |29] l for 
arbitrary densities leads to deviations from the linear density dependence as shown in Fig. [2] In contrast to the 
bound states, there is no Pauli shift for the scattering states. The lowering of the continuum edge is given by the 
Fock energy shift at zero momentum, 

1 /3 

AS^-^p = 0) = - E ni)fe{q) = = -4 (^) . (32) 

q ^ ^ 

The Fock shift at the Fermi momentum AE^°'^^{pp) results as AE^°'^^{p = 0)/2. This is in agreement with 
the estimate of the shift of the continuum by the chemical potentials given above since in the strongly degenerate 
case the chemical potential of the electrons is given by the energy at the Fermi momentum. The lowering of the 
continuum edge is also shown in Fig. [2] The bound states vanish when they meet the continuum edge. We see 
that for n' < 0.0018, only two levels survive and for n' < 0.00018 still 3 levels are left. Only at very small 
densities more than 3 levels are left. We repeat here that the solution of an in-medium Schrodinger more complex 
than Eq.([28ll was dealt with in |[ll[3l[36l[37l|38l[39l[40l[4T]. In |[39l it was shown that perturbation procedures 
are justified even for in-medium Schrodinger equations. Solutions beyond perturbation theory can also be found 
within a variational approach 1161 l42l . 



Fig. 2 The shifts of the ground state and two excited states in Pauli-Fock approximation in relation to the lowering of the 
continuum level. 

As an application we calculated the pressure related to the classical ideal pressure for the case of high densities 
at a temperature of 10000 — 20000 K, corresponding to the Hartree-Fock and Wigner approximations, see Fig. [3] 
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The energy shifts are taken into account. For comparison we have given a curve for 10 000 K given by Vorberger 
et al. P3l and with a resuh for T = 20000 K obtained by means of a Pade approximation in combination 
with a mass action law |19|. The agreement of our new, rather simple approach with the previous ones is rather 
satisfactory. 




23.4 23.6 23.8 24 24.2 24.4 246 24.S 25 23.2 23.4 23.6 23.8 24 242 244 24.6 248 25 

IgditI, T = 10000; 15000; 2000O K Ig(nf), T = 20000 K 

Fig. 3 Calculations of the pressure corresponding to the Hartree-Fock approximation taking into account the energy shifts 
for the case of high densities. The left panel shows p/pid for the temperatures T = 10000, 15000, 20000 K. For logn < 23.5 
begins the region of developed bound states which is not well described by this method. The magenta curve represents 
numerical and analytical results by Vorberger et al. 1431 for T — 10000 K. In the right panel we show the pressure related 
to the classical ideal pressure at the same temperatures and for comparison we give for T — 20000 K a curve (in magenta) 
obtained earlier using a Fade approximation for the pressure in combination with a nonideal Saha equation |19| . 



4 Discussion and Conclusion 

We discuss here the divergence of the partition function of the Bohr atom and the BPL method of renormalization, 
and present several new results within the physical picture. At small densities we discussed closed expressions 
for the ring sums. At high densities we concentrated on the Hartree-Fock approximation and took into account 
the shifts in Hartree-Fock approximation. We assume that these shifts provide the most important effects for the 
destruction of bound states in the high-density region. In particular we contribute here to the theory of hydrogen 
at high pressures in the region where a Mott transition to full ionization has been predicted and where recent 
experiments have shown a transition from insulating behavior to metal-like conductivity. 

In order to understand this transition several effects have to be taken into account. We concentrated here on 
so-called Pauli blocking effects expressing the rule that states occupied by atomic electrons cannot be occupied 
by free electrons with the same spin state. This leads at high electron densities to the destruction of atomic states 
which need a relatively high amount of phase space. We calculated the energy shifts due to Pauli effects and 
discuss the Mott effects solving effective Schrodinger equations for strongly correlated systems. 

Still a word on the plasma phase transition. In the temperature region 10000 — 20000 K which we studied 
here, no first order transitions were observed. The physical picture which we used here cannot be extended 
without problems to lower temperatures, where H2 molecules dominate. In order to check for phase transitions 
we introduced the energy levels calculated here into the chemical picture developed earUer lfT6l which includes 
H2 molecules. A resulting isotherm for T — 6000 K shows a small wiggle. According to this estimate the 
transition occurs in the range below 6000 K, the critical temperature Tc is slightly above. Summarizing, we show 
here that in the non-degenerate the most essential effects are the screening and the BPL reduction, resulting in 
absorbing the orders in e^, and the compensation effects of the states below and above the series limit. Another 
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important effect is the summation of an infinite series in the parameter 7 = A^(TBPh{T) leading to a Saha-type 
behavior At high density the PauH-Fock contributions are most essential for the limitation of the number of 
energy levels and determine the transition from partial ionization to full ionization when densities become very 
high . 
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